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Abstract 

A superconducting rod with a magnetic moment on top develops vortices obtained here through 3D 
calculations of the Ginzburg-Landau theory. The inhomogeneity of the applied field brings new prop- 
erties to the vortex patterns that vary according to the rod thickness. We find that for thin rods (disks) 
the vortex patterns are similar to those obtained in presence of a homogeneous magnetic field instead 
because they consist of giant vortex states. For thick rods novel patterns are obtained as vortices are 
curve lines in space that exit through the lateral surface. 

Key words: Vortex pattern, Ginzburg-Landau, Magnetic dot 
PACS: 74.20.-z, 74.20.De 



p 1. Introduction 



P 

O 

o 



> 
in 

i— H 



o 
o 



1 



The minimal condition for the onset of a vor- 
tex in a mesoscopic superconductor depends on 
geometrical parameters. A thin disk of radius R 
(thickness D <~ £, £ is the coherence length) can 
only exist in the Meissner state for R ~ £. But 
for £ < R < 2£, giant vortex states are allowed 
and for R > 2£ multivortex states become pos- 
sible, as reported in [1] and recently observed in 
Refs. [2, 3]. Mesoscopic superconductors have 
new and interesting properties [4], and also pro- 
vide an interesting playground to understand the 
co-existence of magnetism and superconductivity 
[5, 6]. For instance magnetic dots on top of su- 
perconducting film have been investigated both 
theoretically and experimentally [7-10]. 

In this paper we report a theoretical study done 
on superconducting rods of radius R and varying 
thickness D ) with a magnetic dot on top, a sys- 
tem which displays curved vortices triggered by 
the inhomogeneity of the magnetic field. We ad- 
dress the minimal geometrical conditions for the 
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onset of vortices and also the nature of the vortex 
patterns. 

The shape of a vortex shifts from a flat coin 
to a curved line in three-dimensional space, as 
the thickness D tales the rod from a disk to a 
tall rod. Consequently the vortex pattern con- 
sists of top to bottom giant vortices and top to 
side mult i- vortices in these two extreme limits, 
namely, D ~ £ and D>^, respectively. Rods of 
intermediate thickness display highly non trivial 




Figure 1: (Color online) Superconducting rod with radius 
R and thickness D. The magnetic dot is oriented along 
the x-axis and positioned 2£ above the top surface. 
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vortex patterns showing features of both extreme 
limits. Our results are obtained in the context 
of the Ginzburg-Landau theory in the limit of no 
magnetic shielding. We take a magnetic dot with 
magnetic moment \± positioned 2£ above the rod's 
top surface and oriented along the rod's axis, see 
Fig. 1. 

Although curve vortices also appear in presence 
of a homogeneous field tilted with respect to the 
surface, as vortices must always emerge perpen- 
dicular to the surface [11], the inhomogeneity of 
the field brings new effects. The dipolar magnetic 
field weakens with distance and may become so 
dim at the bottom of the rod to the point that this 
region becomes unable to sustain vortices. Conse- 
quently the Meissner state is kept at the bottom 
of the rod. This favors a reentrant vortex state, 
as a top to side vortex must disappear at the top 
surface, the same one of its onset. Thus the Meiss- 
ner phase is retrieved at a higher magnetic mo- 
ment. This qualifies superconducting rods with a 
magnetic moment on top for technological appli- 
cations as logical switchers, since vortices can be 
expelled or accepted by fine tuning the external 
magnetic moment dot strength. 

2. Theoretical approach 

According to the standard Ginzburg-Landau 
approach, the Gibbs free energy of the super- 
conductor near the critical temperature T c can 
be expanded in powers of the complex order pa- 
rameter \^(r) that gives the density of Cooper- 
pairs: |^(r)| 2 . Hence the Gibbs free energy differ- 
ence between the superconducting and the normal 
states is, 

F s -F n = J dV {a(T)\V\ 2 + + 
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with the phenomenological constants a(T) = 
a (T — T c ) < 0, /3 > 0, m* (the mass of one 
Cooper-pair), and $ = hc/2e, the fundamental 
unit of flux. Boundary conditions are imposed to 
the external surfaces of the rod. We express quan- 
tities in this theory in dimensionless units defined 



by the following reduced units: the coherence 
length f(T) = [-h 2 /2m*a(T)} 1 / 2 (lengths), H c2 = 
0o/27r£(T) 2 (magnetic field), H c2 £ = $ /2<(T) 
(vector potential), /x = H c2 ^{Tf = $ £(T)/27r 
(magnetic moment), = \/ — a (T)//3 (order pa- 
rameter), and F = a{T) 2 /2/3 (free energy). The 
magnetization is calculated using the expression 
M = f fx J(r)dv. 

In terms of these dimensionless units, the Gibbs 
free energy difference becomes, 
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The integration is restricted to the volume of the 
rod. The boundary condition in dimensionless 
units becomes, 



n 



(V - iA)^ 



boundary 

For the magnetic dot we use field produced by a 
point-like dipole A = (jH x f) /r 3 . 

3. Numerical results 

We solved Eqs. 1 and 2 using the simulated an- 
nealing procedure [4]. We set the following pa- 
rameters for the rod: the radius R ranges from 1£ 
to 4£; the thickness D ranges from 2£ to 8£ and 
the magnetic moment \i ranges from to 100//o- 
For the thinnest rod considered here, D = 2£, we 
find that the inhomogeneity of the field does not 
matter, since the free energy and the magnetiza- 
tion are similar to the homogeneous applied field 
case. Clearly this is because the dot's magnetic 
field does not vary significantly inside the rod. 

For R = 1£ and any thickness we find no vor- 
tex state, thus the Meissner state prevails up to 
the normal state. We do not observe ant i- vortices 
or Meissner state with opposite orientation in the 
rods because we limit the present study to rods 
with a maximum radius of 4£. Thus only giant 
vortex states (GVS) are observed. For R = 2£ 
and for D = 2£, only one vortex is allowed. For 
larger R more than one vortex is possible and 
their number increases very rapidly with the ra- 
dius. Although there is more space to accom- 
modate extra vortices the magnetic field becomes 
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Figure 2: (Color online) Free energy F/Fq and magnetiza- 
tion M/Mq versus magnetic moment /i//io f° r mesoscopic 
rods with thickness D = 2£ and several radii. Giant vor- 
tices with increasing maximum vorticity are observed by 
increasing the radius. 



Fig. 3 shows the free energy for the R = 4£ and 
D — 6£ rod. 

The 3D view of the Cooper pair density, |^| 2 , 
shows curved top-to-side vortices, that is, vortices 
that enter from the top surface and exit perpen- 
dicularly in the side surface. We call them N-fold 
states, where N describes the number of vortices. 
Along the energy curve vs /i in Fig. 3 the states 
are labeled according to this definition. The N- 
fold states are not clearly seen in Fig. 3 because 
their free energy lines is superposed to the Meiss- 
ner line from where is their onset and disappear- 
ance. However these states are clearly seen to 
exist by analyzing the \ip\ 2 isosurface, as shown 
in Fig. 4. In the two biggest rods considered here 
we observed the nucleation of N-fold states fol- 
lowed by Meissner and GVS. This features are 
justified by the growth of a normal region closer 
to the the magnetic dot position. When this nor- 
mal region occupies a significant part of the rod 
the remaining superconducting part behaves as a 
thin rod, i.e., it does not exhibits N-fold states. 
So, we retrieve the Meissner state. This feature is 
in Table 1 and correspond to the second Meissner 
state. 



weaker at distances away from the center. Conse- 
quently the average magnetic magnetic field does 
not grow with the radius. This effect is equivalent 
to that studied by Ovchinnikov sometime ago [12]. 
A sawtooth behavior in the magnetization with 
bigger H c2 is clearly established for bigger radius. 
This behavior is observed for R = 3£ and D = 2£. 
In Fig. 2 we show the vortex states obtained for 
different thin rods, i.e., disks. Their magnetiza- 
tion curves are similar to those obtained in pres- 
ence of homogeneous field. In table 1 we give the 
complete description of the states through the set 
of parameters (i?, D, /i//x ). 

The inhomogeneity of the field becomes impor- 
tant for the thicker rods D = 6£ and 8£, as shown 
here. Even for small radii, the rod already dis- 
plays behavior that differs significantly from the 
homogeneous field case. We find smooth free en- 
ergy lines, with no first order transition occa- 
sioned by the entrance of quantized vortex lines. 



R=4^ and D=8^ 
o,o -| 




Figure 3: (Color online) Free energy F/Fq vs fi/fio for a 
rod with R = 4£ and D = 8£. For increasing magnetic 
moment the vortex pattern evolves continuously through 
the N-fold configurations, as labeled in the inset. See, 
Fig. 4 for the corresponding |^| 2 isosurface views of these 
N-fold states. 
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Table 1: The sequence of ground states for the of rods considered here. The magnetic dot is positioned 2£ away from 



the rod's top surface and oriented along the z-axis. The question mark means the normal state was not reach. 



R/Z 

D/Z 
State 
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1 
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1 
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2 
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2 
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2 
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2 
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3 
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3 
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3 
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3 
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4 
2 


4 
4 


4 
6 


4 
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Meissner 


0-36 


< 100 


< 100 


< 100 


0-18 


0-76 


<100 


<100 


0-14 


0-46 


<100 


<100 


0-14 


0-20 


0-22 


0-22 


1 GVS 










20-30 


78-100 






16-24 


48-76 






16-22 


22-100 






2 GVS 










32 








26-30 


78-100 






24-30 








3 GVS 


















32-38 








32-36 








4 GVS 


















40-50 








38-44 








5 GVS 
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3-fold 




























34-38 


24-34 
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4-fold 




























40-44 


36-42 


34-40 


5-fold 
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6-fold 
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50-58 


7-fold 
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Meissner 






























64-74 


60-100 


1 GVS 






























76-100 




Normal 


> 36 
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> 32 
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> 50 
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? 
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? 



4. Conclusions 

Using Simulated Annealing, a truly three- 
dimensional approach, we obtain the vortex pat- 
terns of mesoscopic rods with a point-like mag- 
netic moment on top. We conclude that rods with 
thickness smaller than 4£ can be considered as 
thin disks since only top to bottom giant vortex 
states are obtained. For rods with thickness big- 
ger than 4£ , we observed the appearance of N-fold 
vortex states, namely N top-to-side vortices. The 
onset and disappearance of these states is from 
a single GVS or Meissner state line, as contin- 
uous transitions from it. Consequently, for suffi- 
cient large magnetic moments one can retrieve the 
Meissner state before reaching the normal state, 
and this, leads to a reentrant behavior. 
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